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Abstract 

A covariant Lagrangian formulation of a solution to the cosmological constant 
problem, based on vizualising the fluctuations of the vacuum energy as a non- 
equilibrium process with stochastic behaviour, is presented. The variational 
principle yields equations of motion for the cosmological "constant" A, treated 
as a dynamical field, together with an equation for a Lagrange multiplier field 
4>, and the standard Einstein field equations with a variable cosmological 
constant term. A stochastic model of A yields a natural explanation for the 
smallness or zero value of the constant in the present epoch and its large value 
in an era of inflation in the early universe. 
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A recent model for solving the cosmological constant problem has been proposed in 
which the vacuum energy is treated as a fluctuating environment with stochastic behaviour. 
In the following, we shall present a covariant formulation based on a Lagrangian density, 
which yields classical equations incorporating Einstein's gravitational field equations, upon 
which a stochastic treatment using a Wiener process can be developed. 

Although the model uses methods of critical phenomena and non-equilibrium statistical 
mechanics to model the vacuum energy, it can be considered as a phenomenological de- 
scription of the kinds of behaviour that could be expected in a more fundamental quantum 
gravity theory. It is generally agreed that the cosmological constant problem cannot be 
solved within the context of a purely classical theory of gravity. However, no satisfactory 
quantum gravity theory has been formulated, so it is hoped that our model can shed light 
on the solution to the problem without the full apparatus of such a theory. 

The Lagrangian density is given by 

C = C R + C A + C M , (1) 

where 

£r = yf^sTR^ (2a) 

C A = -2v^[A + (A,X - a A + A 2 )0] , (2b) 

and A = A(x) is the variable cosmological "constant", treated as a dynamical field, is a 
Lagrange multiplier field, m m = dx^/dr is an observer's four- velocity along a world line in 
spacetime, a is a constant and Cm is a matter Lagrangian density. A variation of C with 
respect to <fi and A yields the equations of motion: 

A iA y - aA + A 2 = 0, (3a) 
y=^V>), M + (a-2A)0-l = O. (3b) 



Varying C with respect to g^ u and using (|3"a]) gives 



R^u ~ ^9^uR + Agr^ = 8-itGT^. (4) 



The kinematical variable is the four-velocity of a fluid element along a world line in 
spacetime, associated with a fluid with density p and pressure p, so we do not vary C\ with 
respect to u/\ 

The cosmological constant enters through the vacuum energy density: 

Tvpv = —Pvdnv = ~~^q9^v- (^) 

Today, A has the small value, A < 1CT 46 GeV 4 , whereas generic inflation models require 
that A has a relatively large value during the inflationary epoch. This is the source of the 
cosmological constant problem. 

The line element in the Friedmann-Robertson- Walker model is 

dr 2 = dt 2 - R 2 (t)\^— + r 2 d6 2 + r 2 sin 2 9d<f> 2 }, (6) 
1 1 — kr 2 ) 

where k = — 1,0, +1 and we have used comoving coordinates with = (0,0,0,1). Then, 
Eqs.Q3aD, © and (§ become 

A = aA - A 2 , (7a) 

0+— 0+(a-2A)0-l = 0, (7b) 

where A = dA/dt, pu denotes the mass density, H is the Hubble constant, and in the 
following H and to denote the present values of H and t, respectively. We define 

^tot = ^ M + = 1 — fife, (8) 

where fi = 8nGp/3H 2 . 

We shall treat the vacuum energy as a fluctuating environment and consider A as a 
variable characterizing the state of this system. The parameter a in Eq . (|7aD corresponds 
to the difference between the growth and decline of particle-antiparticle annihilation in the 
vacuum, while the second term is a self-restriction term which limits the growth of A. 
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In ref.(l), we considered the situation in which the vacuum fluctuations are rapid com- 
pared with Tmacro = a ~\ which defines the macroscopic scale of time evolution. We 
assumed that the parameter a can be written as at = a + cr£t, in which a is the average 
value, £t is Gaussian noise and a measures the intensity of the vacuum fluctuations. Let us 
write Eq.(|7ej) as 

dAt = (aA t - A 2 t )dt + aA t dW t = f(A t ) + ag{A t )dW u (9) 

where dW t is a Wiener process. The probability density p(x, t), which describes the A 
distribution, satisfies the Fokker-Planck equation: 

d t p(x,t) = -d x [(ax - x 2 )p(x,t)] + —d xx (x 2 p(x,t)). (10) 

The diffusion process is restricted to the positive real half line and and oo are intrinsic 
boundaries, because g(0) = and /(oo) = — oo. The probability of the diffusion process 
reaching infinity as t — > oo is zero, since infinity is a natural boundary. Moreover, zero is a 
natural boundary if a > cr 2 /2, so neither boundary is accessible and no boundary conditions 
need be imposed on the Fokker-Planck equation. For a < cr 2 /2, it can be shown that zero 
is an attracting boundary. 

The stationary-state solution for the probability density, p s (x), of Eq. ([To|) is given by 

p s (x) = Nx^^- 2 ex P (~^). (11) 



The normalization constant iV is 



N- 1 



2 v2(a/«r»)-l-,-l (2a 



a 2 



r 



a 2 



where T denotes the T-function. If p(x,t) is integrable between and oo, then a stationary 
state solution exists when a > a 2 jl. If it does not exist, then the probability density will 
be concentrated at zero, i.e., p(A) = 5(A) for a < a 2 /2. 

For < a < a 2 /2, the vacuum fluctuations dominate over the growth or decline of A, 
although the value zero is still the most probable value for A, since the distribution function 



has a vertical slope at A = 0. Because we are using a continuous variable, A never reaches 
the boundary zero in a finite time. 

When a > cr 2 /2, the growth of A dominates the influence of the vacuum fluctuations, 
and in the neighborhood of zero the probability of A = drops to zero. For the stochastic 
model there are two transition points described by different order parameters. At a = cr 2 /2 
real growth of A becomes possible corresponding to a change from a degenerate random 
variable for steady-state behavior to a stochastic variable; the boundary at a = switches 
from attracting to natural. Secondly, there is the transition point a = a 1 which corresponds 
to a qualitative change in the stochastic variable A with no change in the nature of the 
boundary. The probability of A = drops abruptly to zero. 

The following scenario can be deduced from our model. In the inflation era, the intensity 
of vacuum fluctuations is large and a > a 2 , causing a second-order phase transition and a 
maximum in A not near zero. This corresponds to the large vacuum energy needed to drive 
inflation || . As the universe expands the intensity of vacuum fluctuations decreases and for 
< a < cr 2 /2 or er 2 /2 < a < a 2 the probability density is largest when A is non- vanishing 
and small, which can lead to a current value of A that can be used to fit the observational 
data. If the stationary probability density p s does not exist for a < <J 2 /2, then A = is a 
stationary point; the drift and diffusion vanish simultaneously for A = and p(A) = 6(A). 
This corresponds to the case when A is vanishingly small. 

Thus, our model provides a natural explanation, in terms of non-equilibrium stochastic 
processes in an expanding universe, for the behavior of A required to fit observational data 
and still be consistent with inflationary models. 

Ongoing searches |4[] for Type la supernovae show that ^a < 0.47 (at 95% confidence 
level for spatially flat A models). Moreover, for A ^ a larger fraction of QSOs would be 
gravitationally lensed and QSO surveys give Q\ < 0.7 |J. Cold dark matter models (CDM) 
for large scale stucture formation, which include a cosmo logical constant, yield a better fit 
to the shape of the observed spectrum of galaxy clustering than does the standard Qm = 1 
CDM models, using h = i/o/(100km/sec/Mpc) = 0.7, Q\ = 0.6 and a baryon density with 
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Qb = 0.0255, consistent with primordial nucleosynthesis ||. However, the amplitude for the 
A CDM models is too high compared to the data, a problem that persists at all scales. 

The problem of the age of the universe is also alleviated in A models. An analysis of 
the cosmological data showed that for Q\ = 0.65 ± 0.1, Qm = 1 — and a small tilt: 
0.8 < n < 1.2, models exist which are consistent with the available data and an inflationary 
spatially flat universe J?]]. 

Models based on A treated as a scalar dynamical field || have been found to partially 
resolve observational problems. The observations of gravitationally lensed QSOs yield a less 
restrictive upper bound on Hot® in these models They may also provide a solution to 
the size of the amplitude problem, since although the shape of the spectrum is the same as 
that of the A CDM model with = 0.6, the dynamical A model yields a lower amplitude 
and therefore gives a better fit to the galaxy clustering data || . 
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